Void-Supercluster Alignments 

Daeseong Park and Jounghun Lee 

Department of Physics and Astronomy, FPRD, Seoul National University, Seoul 151-747, 

Korea 

pds2001@astro . snu . ac . kr , j ounghun@astro . snu . ac . kr 

ABSTRACT 

We investigate alignments between the spin axes of cosmic voids and the 
principal axes of nearby superclusters using the Millennium Run simulation of 
a ACDM cosmology. The concept of void spin was first introduced by Lee and 
Park in 2006 to quantify the tidal effect on voids from the surrounding matter 
distribution. Our numerical analysis reveals that the void spin axes are strongly 
aligned with the supercluster minor axes, but anti-aligned with the major axes, 
and have no correlations with the intermediate axes. We provide physical expla- 
nations to this numerical results on the basis of tidally induced correlations. It 
is expected that our work will provide a new insight into the characterization of 
the cosmic web on the largest scale. 

Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 



Voids stand for those regions in the universe which have extremely low density. There 
is still no consensus on how to define voids. In recent numerical approaches using N-body 
simulations, voids are often defined as regions empty of massive dark matter halos. In 
real observations, however, voids are usually defined as und erdense regions with the mean 



density contrast close to —0.9 (e.g., iHoyle fc Vogeleyl 12004 ). given the limitation that we 



cannot observe very faint galaxies. Adopting the latter definition of voids, one can say that 
voids are not empty of massive dark halos but have some structures in them. An interesting 
question is how the void structures evolve with time. 



The standard picture based on the cold dark matter paradigm explains that the voids 
originated from the local minima of the initial density fluctuation and have become more and 
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more underdense via gravitational rarefaction in the subsequent evolution (Ilckdll984j ). Re- 
cent high-resolution simulations, however, have reveale d that the evolution of the void struc- 
tures is more complicated than naively predicted (e.g., iGottlober et al.ll2003l ; IColberg et al. 



20051 ): Not all voids remain underdense but some fraction of them actually collapse into 
bound halos. Although the collapse of voids may be described qualitat ively as the occurrence 
of the clouds-in- voids in the frame of the standard excursion set theory (ISheth fc van de Weygaert 
2004 ). it has yet to be well understood under what conditions and how frequently the collapse 
of voids takes place. 

A first suspic ion goes to t he tidal forces f r om the sur rounding matter distribution 
(jSahni et al.lll994j ). As noted by lShandarin et al.l (12004 120061 ). the voids should be severely 
disturbed by the tidal influence due to the low density. Considering this possibility, a plau- 
sible scenario is that when the tidal effect on a void wins over the gravitational rarefaction, 
the void will become squeezed out to collapse into a bound halo. A remaining question is 
how to quantify the tidal effect on voids and predict its consequences. 



Very recently, iLee fc Park! (120061 . hereafter, LP06) introduced a new concept, the void 
spin angular momentum, which is defined as 

Nv 

J = -rr- ^m Q r Q X v Q , (1) 



1 

My 



where N v is the total number of halos belonging to a void, m a , r a , and v Q are the mass, the 
position, and the velocity of the a-th halo belonging to a void and My = J2a V m c 



Here 



the positions and the velocities of the void halos are measured relative to the center of mass 
of the void halos. 

As mentioned by LP06, although the quantity defined in equation (TjQ) is not a real spin 
angular momentum since voids are not bound system, the concept of void spin turned out 
to be a very effective measure of the tidal effect from the surrounding matter. To show that 
the void spin angular momentum is tidally induced, LP06 measured the spatial correlations 
between the spin axes of neighboring voids analyzing numerical data from high-resolution 
simulation, and found that the numerical results are in excellent agreement with the analytic 
predictions based on tidally generated correlations. 

If the void spin angular momentum is really induced by the tidal effect from the sur- 
rounding matter, then there should exist correlations not only between the spin axes of 
neighboring voids but also between the spin axes of voids and the neighboring matter dis- 
tribution. This cross correlation, if found to exist, should be a more direct indication of the 
tidal effect on voids from the surrounding matter. 

Our goal here is to measure the cross-correlations between the void-spin axes and the 
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axes of the neighboring superclusters using numerical data from high-resolution simulations 
and to explain it physically. The superclusters are considered as the counterparts of voids 
since they are the largest bound structures, comparable to voids in size. 

The outline of this paper is summarized as follows. In §2, the numerical data are 
analyzed and the signals of void-supercluster alignments are detected. In §3, analytic model 
for the void-supercluster alignments is derived and compared with the numerical results. In 
§4, the implications and caveats of our work are discussed, and a final conclusion is drawn. 



2. NUMERICAL RESULTS 
2.1. Measurements of Void Spins 



We analyzed the halo catalog extracted from the Millennium Run simulation (ISpringel et al 



20051 ). which is now available at http://www.mpa-garching.mpg.de/millennium, The sim- 
ulation is of a flat ACDM universe for which the values of the key parameters are given 
as (fl m ,a 8 ,h,n s ) = (0.25,0.9,0.73,1) where fl m , cr 8 , h and n s represent the matter density 
parameter, the rms density fluctuation on the top-hat scale of 8/i -1 Mpc, the dimensionless 
Hubble parameter and the slope of the primordial power spectrum. 

For our analysis, we exclude those halos from the Millennium Run catalog which have too 
low particle numbers, as those low-particle number halos are just poorly sampled, heavily 
affected by noise (V. Springel in private communication). Technically, we set the cut-off 
particle number at 50, a minimum particle number for defining a halo density profile. 



In accordance with the void-finder algorithm proposed by iHoyle fc Vogeleyl (120021 . here- 
after, HV02) to the Millennium Run halo catalog, we take the following steps to identify 
voids: First, we classify the halos according to the wall/field criterion, given as / = d 3 + 3a/2, 
where d% is the average distance to the third nearest neighbor and a is its standard deviation. 
For the Millennium Run halos, the criterion distance is found to be Z = 2.78/i _1 Mpc. Note 
that this value is s lightly higher than the value I = 2.44/i _1 Mpc used in our previous work 



(ILee fc Park! 120061 ) . It is because in our previous work the voids were identified not from 
the halo catalog but from the galaxy catalog. Second, we place the wall halos on cubic grid 
cells each of which has a linear size of I, counting the number of halos in each cell. Third, 
increasing the radii of empty spheres from the center of all empty grid cells till each of them 
include three wall halos on the surface, we determine the largest possible empty spheres. 



Finally, we detect voids by categorizing the overlapping empty spheres whose radii 
are greater than the minimum size threshold. The minimum void size are set at 6/i -1 Mpc 
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which was found through statistical significance test in our previous work (ILee h Parkll2006l ). 
A total 20291 voids are identified in the z=0 catalog, among which only 6430 voids are 
found to have more than 30 halos. The mean density contrast (S v ) and the mean effective 
Eulerian radius (Re) of these 6430 voids are found to be 5 V = —0.9 and Re = 12.22/i -1 Mpc, 
respectively. 

Considering only these large voids with more than 30 halos, we calculate the spin angular 
momentum by equation (fTT). It is worth mentioning here that the direction of the void spin 
angular momentum might depend on the cut-off particle number for the halo selection. In 
our analysis, we include only those halos in the Millennium Run halo catalog which contain 
more than 50 particles. One may suspect that a different cut-off particle number might yield 
different directions of the void spins. 

To examine how the void spin angular momentum depends on the cut-off particle num- 
ber, we calculate a second angular momentum J' for each of the 6430 voids, taking all halos 
inside the void with more than 30 particles, and calculate the cosines of the angles, if), 
between J and J' as 

I J • J' I 

cos^ = TWjTr ( 2 ) 

Figured] plots the probability distribution of cosip as histogram with the Poisson errors. The 
horizontal line corresponds to the case of no correlation. As can be seen, the two angular 
momentum vectors are very strongly aligned. It indicates that excluding those halos with 
lower particle number will not affect the void spin angular momentum significantly. 



2.2. Measurements of Supercluster Principal Axes 



A next task is to ident ify superclusters f rom the Millennium Run halo catalog. Following 
the common method (e.g.. lWray et al.ll2006l ). we define a supercluster as a cluster of clusters, 
and identify superclusters with the help of the friends-of-friends algorithm (FOF). The cluster 
halos are se l ected as those halos whose mass exceeds a typical poor cluster mass, M c . As in 
Wrav et all J2006h . we set the value of M c at 1.75 x 1O 13 /i" 1 M . 



Before applying the FOF algorithm to the Millennium Run halo catalog, one has to 
determine a linking length for the supercluster identification. To find an optimal linking 
length, L, we perform a statistical significance test by generate 100 random Poisson samples 
of clusters which have the same number density in the same the box size as the Millennium 
halo catalog but having no clustering effect. The statistical significance of a supercluster can 
be estimated as P(L) = 1 — N po (L)/N sc (L), where N po (L) and N SC (L) are the numbers of 
superclusters found in the random samples and in the Millennium Run simulation data at a 
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linking length L, respectively. (IBasilakos et al.l 120061 ) . 



Figure [2] plots the statistical significance P(L), which reveals that around L = 6/i -1 Mpc 
the statistical significance of finding a supercluster reaches 99% confidence level. Therefore, 
we set the linking length for the FOF algorithm at L = 6/i -1 Mpc which corresponds to the 
linkage parameter b = 0.35. 

A total of 4014 superclusters are found. Among them, a total of 382 superclusters are 
found to consist of more than 5 members. The mean mass of these 382 superclusters are 
found to be M s = 4.2 x 10 14 /i _1 M Q . Using these superclusters, we measured the inertia 
momentum tensor of each supercluster. Rotating the system into the principal axis frame, 
we find the three eigenvectors of the inertia momentum tensor of each supercluster. 



2.3. Alignments between Void Spin and Supercluster Principal Axes 

Now that the spin axes of voids and the principal axes of the superclusters are all 
determined, we calculate the squares of the cosines of the angles between the the void spin 
axes and the three axes of the superclusters as a function of the separation distance, r. Let 
y°) y j y 7 be the three unit eigenvectors of the supercluster inertia momentum tensor, which 
represent the supercluster major, intermediate, and minor axes, respectively. Basically, we 
calculate the following three correlation functions using the selected void-supercluster pairs: 

W «(r) = (|J-yT>(r)-i (3) 
u/(r) = (|J-yT>(r)-^ (4) 
^(r) es (\J-r\ 2 )(r)-l- (5) 

If there is no alignment, then these cross-correlation function will be just zero. 

To examine how these void-supercluster cross correlations depend on the linking length, 
we repeat the whole process using different values of L. The three cross-correlations functions 
between the void spin axes and the supercluster major (uj a ), intermediate (or), and minor 
axes (a; 7 ) for the fives different cases of the linking lengths are plotted in Figs. [31 0] and El As 
can be seen, for those values of L that corresponds to higher than 90% confidence levels (i.e., 
L < 7/i _1 Mpc), there are clear consistent anti-alignment and alignment signals between 
the void spin axes and the supercluster major and the minor axes, respectively, within a 
distance of 30/i -1 Mpc, and there is consistently no alignment signal with the supercluster 
intermediate axes. 
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We provide physic al explanations to this phenomena in §3 within the analytic framework 
proposed originally by lLee fc Pen! (120011 ). 



ANALYTIC PREDICTION 



3.1. Review of Spin-Shear and Direction-Shear Correlations 



According to the Lee- Pen formalism based on the linear tidal torque theory (IDoroshkevich 



19701 ; IWhitelll984j ). the unit spin vector J = ( jj) of a bound halo is correlated with the unit 



traceless tidal shear tensor T = (TL-) as 



{jJj\T) 



1 + a 



5ij - af ik (x)f kj (x) 



(6) 



where a is the spin-shear correlation parameter in the range of [0,3/5]. It represents the 
strength of the correlation between J and T: If a = 3/5, the correlation is strongest. If 
a = 0, there is no correlation. Here, the unit tidal tensor T is intrinsic and local, defined at 
the halo position, x. 



Lee fc Pen! ( 120011 ) also suggested the following formula for the correlation between the 
direction vector to the nearest neighbor y = (y$) with the unit traceless tensor T as 



(M?I T > 



1 



-5ij + 6T ife (x)T fci (x) 



(7) 



Here the parameter, b in the range of [—1, 1] represents the strength of the correlation between 
y and T. Note the difference in the range between the two correlation parameters, a and b. 
The maximum value of b is 1 while that of a is 3/5, less than unity. It reflects the fact that 
a perfect alignment is allowed between y and T but not between J and T. 

Anyway, a crucial implication of equations and (|7D is that the spin axes of halos 
should be closely correlated with the directional geometry of the nearby halo distribution 
since J and y are both correlated with the tidal field. In §3.2, we extrapolate the validity of 
equations (jSJ) and ([7j) which hold good for halos to the voids and superclusters. 



3.2. Modeling Void-Superchister Alignments 



LP06 have already extrapolated the validity of equation (jBj) to unbound voids, assuming 
that the spin-shear correlation parameter a has the maximum value of 3/5 for the case of 
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voids whose spin is denned as ([TJ): 

8 3 

(JiJj\T) = —5ij - -f ik (x v )f kj (x v ), (8) 

where T is now defined at the void center, x„. 

Now, we attempt to extrapolate the validity of equation (J7J) to the alignments between 
the void spin axes and the supercluster principal axes. The superclusters are conspicuously 
elongated along local filaments where the dark matter are preferentially located. The nearest 
neighbors are most likely to be found in the direction along the supercluster major axes. But, 
the filaments are one dimensional structure, collapsed along the major and the intermediate 
principal axes of the local tidal tensors. Thus, the direction of local filaments (i.e., the 
supercluster major axes) are in fact anti-aligned with the tidal tensor major axes, and aligned 
with the tidal tensor minor axes. 

Using the above logic, we assume the following: 



For the direction of the supercluster major axis, y a , the direction-shear correlation 
parameter b has the minimum value of —1: 



(y?y?\T) = -8^ -f ik (x s )f kj (x s ). 



(9) 



For the direction of the supercluster intermediate axis,y^, the direction-shear correla- 
tion parameter b has the value of 0: 



(10) 



• For the direction of the supercluster minor axis, y 7 , the parameter b has the maximum 
value of 1: 

(y]y]\T)=f lk (x s )f kj (x s ), (11) 

Note that in equations (|9l)-(|TTT). the unit tidal tensor is defined at the supercluster center, 
x s . 

Using equations (!B])- (jTT|) . one can derive analytically the three correlation functions, 
u> a , ujP and u 7 , defined in §2.. In this deriva tion, the key part is to calculate the four point 
shear correlation, (T(x 1 ,)T(x 1 ,)T(x s )T(x s )). iLee fc Pen! ( 20011 1 calculated this quantity for 
the case that the two unit tidal tensors are defined at the same halo position but smoothed 
on two different scales. What they found is t hat it is approximated as the density auto- 
correlation (see Appendix I in lLee &: Peril l200ll ). In our case, the two unit tidal tensors are 
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not defined at the same position but at two different positions, the void and the supercluster 
centers. Therefore, the four-point shear correlation would not be approximated as the density 
autocorrelation. Instead, it may be approxi mated as the dens ity two-point correlation, just 
like the void spin-spin correlation function (ILee fc Parkl 120061 ) . 



Hence, using the same approximation used in (ILee fc Parkl|2006l ) but considering the 
fact that the maximum value of the correlation parameter b is different from that of a, we 
find that the void-supercluster alignments can be approximated as 



uo [r) 



i e 



R\ 



io £i(o) 



1 &(r) 

io£i(o) : 



(12) 
(13) 
(14) 



where r = |x„ — x s |. Here is the two point correlation function of the density field on 
the smoothing scale of R. For the void-supercluster alignments, the smoothing scale should 
be a minimum Lagrangian radius enclosing a void-supercluster pair. Since the separation 
distance r cannot decrease below the sum of the supercluster radius and the void diameter 
in Lagrangian space, the smoothing scale may be written as 

R = R s + 2R V , (15) 

where R s and R v represent the effective radii of a supercluster and a void, respectively. 

For the comparison with the numerical results, we relate the Lagrangian supercluster 
and void radii to the observables given in §2: 



R« 



1/3 



^ = (1 + 4) 1/3 ^, 



(16) 



where p is the mean mass density of the universe. 



Equations (fT2"|) . f[T3"j) and (TH|) are plotted in Figs. El CD and El respectively, where the 
numerical results (with L = 6/i~ 1 Mpc) derived in §2.3 are also plotted as solid dots with 
Poissonian errors. For the analytic resu lts, we use the transfer function of the initial power 
spectrum given by iBardeen et al.l (119861 ) with the cosmological parameters set at the values 
used in the Millennium Run simulations. And, the shape parameter, T is approximated as 
T = Q m h (V. Springel in private communication). As can been obviously seen, the analytic 
approximations work very well. The excellent agreements between the analytical and the 
numerical results imply that the void-supercluster alignments are indeed generated by the 
tidal interactions between the voids and the superclusters. 
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4. DISCUSSION AND CONCLUSION 

We have investigated correlations in spatial orientations between voids and their neigh- 
bor superclusters using the Millennium Run simulation of a concordance cosmology. Adopt- 



ing the concept of void spin proposed by lLee &: Park! (120061 ). we have found for the first time 



that the void spin axes are very strongly correlated with the supercluster minor axes within 
the separation distance of 30/i _1 Mpc. Testing how the result depend on the choice of void 
and supercluster definition, we have found that our result is quite solid. 

To this numerical phenomena has a physical explanation been provided based on tidally 
generated correlations. Under the assumption that the neighboring superclusters are rep- 
resentative of the filamentary surrounding matter which wrap and exert tidal forces on the 
voids, we have derived an analytic formula for the alignments between the void spin axes 
and the supercluster minor axes. The analytic prediction has turned out to agree with the 
numerical result very well. 

However, it is worth discussing a caveat which the success of our work is subject to. This 
caveat lies in the fact that there is no consensus on how to define voids and superclusters, 
unlike the case of halo-defining. Using different algorithms could produce different results for 
correlations between voids and superclusters. Although we have shown here that the void- 
supercluster alignments do not strongly depend on the linking length of the FOF algorithm 
and the cut-off mass of the void halos, it will definitely necessary to compare between different 
methods for the void and supercluster identifications in the future. 



Together with our previous work on the void spin-spin correlation flLee fc Parkl 120061 ) . 
this new result supports the scenario that the voids originate from the initial regions where 
the tidal effect becomes maximum. In addition, this new result on the void-supercluster 
alignments demonstrate more directly how the largest scale structure and voids are connected 
in a cosmic web through tidal influences. 

A final conclusion is that our work will provide a new clue to describing the void- 
supercluster network in the cosmic web and leads us to have a deeper insight into the 
formation and evolution of the large scale structure of the universe. 



The Millennium Run simulation used in this paper was carried out by the Virgo Super- 
computing Consortium at the Computing Centre of the Max-Planck Society in Garching. 
We thank V. Springel and G. Lemson for plenty of helps. This work is supported by the 
research grant No. R01-2005-000-10610-0 from the Basic Research Program of the Korea 
Science and Engineering Foundation. 



REFERENCES 

Bardeen, J.M., Bond, J.R.,Kaiser, N., & Szalay, A.S. 1986, ApJ, 304, 15 

Basilakos, S., Plionis, M., Yepes, G., Gottlober, S., Turchaninov, V. 2005, MNRAS, 365, 539 

Colberg J. M., Sheth R. K., Diaferio A., Gao L., & Yoshida N. 2005, MNRAS, 360, 216 

Doroshkevich, A. G. 1970, astrofizika, 6, 581 

El-Ad, H., & Piran, T. 1997, ApJ, 491, 421 

Gottldber, S., Lokas, E. L., Klypin, A., & Hoffman, Y. 2003, MNRAS, 344, 715 

Hoyle, F., & Vogeley, M. S. 2002, ApJ, 566, 641 

Hoyle, F., & Vogeley, M. S. 2004, ApJ, 607, 751 

Icke V. 1984, MNRAS, 206, 

Lee, J. & Park, D. 2006, ApJ, 652, 1 

Park, D. & Lee, J. 2007, submitted to Phys. Rev. Lett. 

Lee, J. & Pen, U. L. 2001, ApJ, 555, 106 

Sahni, V., Sathyaprakah, B. S., & Shandarin, S. F. 1994, ApJ, 431, 20 

Shandarin, S. F., Sheth, J., & Sahni, V. 2004, MNRAS, 353, 517 

Shandarin, S., Feldman, H. A., Heitmann, K., & Habib, S. 2006, MNRAS, 367, 1629 

Sheth, R. K., & van de Weygaert, R. 2004, MNRAS, 350, 517 

Springel, V. et al. 2005, Nature, 435, 629 

Wray, J. J., Bahcall, N. A., Bode, P., Boettiger, C, Hopkins, P. F. 2006, ApJ, 652, 907 
White, S. D. M. 1984, ApJ, 286, 38 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 11 - 



1.5 - 



in 
O 
o 



1 



0.5 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



rr 



J I I I I I I I I I I I I I I I I I L 



0.2 0.4 0.6 0.8 

cos^ 



1 



Fig. 1. — The probability density distribution of the cosines of the angles between J and J'. 
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Fig. 2. — The statistical significance of the linking length for the FOF algorithm to find 
superclusters in the Millennium Run halo catalog. 
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Fig. 3. — Numerical result of the cross-correlation between the void spin axes and the major 
axes of the neighboring superclusters as a function of separation distance for the five different 
cases of the linking length, L. 
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Fig. 4. — Same as Figure [6] but with the supercluster intermediate axes. 
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Fig. 5. — Same as Figure [6] but with the supercluster minor axes. 
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Fig. 6. — Comparison of the analytic cross-correlation between the void spin axes and the 
major axes of the neighboring superclusters (solid line) with the numerical result from the 
Millennium Rum simulation (dots). The linking length for the numerical result is set at 
6 _1 Mpc which corresponds to the 99% confidence level for the supercluster identification 
using the FOF algorithm. The errors represent the standard deviation for the case of no 
alignment. The horizontal line corresponds to the case of no alignment 
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Fig. 7. — Same as Figure [6] but with the supercluster intermediate axes. 
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Fig. 8. — Same as Figure [6] but with the supercluster minor axes. 



